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The one-dimensional 0(2) model is the simplest example 
of a system with topological textures. The model exhibits 
anomalous ordering dynamics due to the appearance of two 
characteristic length scales: the phase coherence length, L ~ 
t 1 ^ , and the phase winding length, L w ~ L x . We derive the 
scaling law z = 2 + fi\, where fi — (/i = 2) for nonconserved 
(conserved) dynamics and x = 1/2 for uncorrelated initial 
orientations. From hard-spin equations of motion, we consider 
the evolution of the topological defect density and recover a 
simple scaling description. 
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The scaling hypothesis has played an important role in 
our understanding of the late-stage ordering dynamics of 
systems quenched from a homogeneous disordered phase 
into an ordered phase region with a broken symmetry 
[1]. According to this hypothesis, the order parameter 
morphology at late times after the quench is statistically 
independent of time if all lengths are rescaled by a sin- 
gle characteristic length scale L(t). This implies that 
the pair correlation function C (r, t) of the order param- 
eter should depend on its arguments only through the 
ratio r/L(t). Recently, we have shown [2] that a natural 
extension of the scaling hypothesis to two-time correla- 
tions, C(r,t,t') — f(r/L(t),r/L(t')), supplemented by 
an understanding of the short-distance (i.e. r <C L(t)) 
structure that follows from any singular topological de- 
fects seeded by the quench [3], determines the late-time 
growth-law of L(t). Any departure from these growth 
laws implies a breakdown of the single-length scaling. 

Given the importance of the scaling phenomenology, 
it is important to look for exceptions to single-length 
scaling, and to try to understand them within a broader 
scaling framework. There is evidence that some systems 
with nonsingular topological textures may violate con- 
ventional single-length scaling. Textures have a spatial 
extent, which can in principle introduce a new character- 
istic length scale. The 0(n) model for an n-component 
vector field in spatial dimension d = n—1 provides a class 
of models with topological textures. Indeed, the 0(3) 
model in d = 2 seems to have at least one of its charac- 
teristic scales growing as t 1 ^ 3 for nonconserved dynamics 
[4,5], which contrasts with the anticipated t 1 / 2 growth 
for this system if scaling holds [2,5]. 

The 0(2) model (or 'XY model') in spatial dimen- 
sion d = 1 is the simplest system with topological tex- 
tures. We show through the time-derivative correlations, 



T(r,t) = dtdf\ t=t , C(r,t,t'), that single-length scaling 
fails in this system, both for conserved and non-conserved 
dynamics. The non-conserved case is exactly soluble 
[6,7]: one finds that C(r, t) scales with a length scale i 1 / 4 , 
different from the t 1 / 2 scaling predicted from the scaling 
hypothesis [2], and T(r,t) scales with the same length- 
scale but with an anomalous time-dependent prefactor. 
We show, within a general framework, that this discrep- 
ancy is due to the existence of two characteristic length 
scales, the 'phase coherence length' L ~ i 1 / 2 , and the 
'phase winding length', L w ~ t 1 / 4 . Since L w is the typ- 
ical length scale over which the phase changes by order 
unity (see Figure 1), it provides the characteristic scale 
for the pair correlation function. The phase coherence 
length, L, drives the dynamics and enters into the two- 
time correlation functions, such as T(r,t). 

In contrast to the nonconserved model, the conserved 
model has not previously been addressed except by com- 
puter simulations. The same concepts are relevant to 
this case also, however, and by means of a simple scaling 
argument we find L ~ i 1 / 3 and L w ~ i 1 / 6 , the latter in 
agreement with recent simulations [7,8]. 

The issue of scaling in these systems is clarified by 
showing that the phase-difference correlation function 
G(r,t) (10), rather than the order parameter correlation 
function, exhibits a generalized form of single-length scal- 
ing with characteristic length L. The topological charge 
density, proportional to the phase gradient, provides an 
equivalent description that exemplifies the simple scaling 
description. 

The key elements of this paper are i) the development 
of hard-spin equations of motion for conserved dynamics, 

ii) the simplification of the equations of motion at late 
times using the separation of the length-scales L and L Wl 

iii) the scaling relations for the length-scales, and the 
consequent forms for spin-spin correlations C(r, t) and 
T(r,t), and iv) the unifying scaling description in terms 
of the topological charge density. 

Figure 1 shows a typical configuration generated by a 
computer simulation with nonconserved dynamics. The 
configuration consists of sections of typical length L 
where the order parameter winds in a given sense, al- 
ternating with antiwinding sections. The winding length 
L w is the typical distance between successive windings of 
2tt in the phase. Since each complete winding (antiwind- 
ing) represents a topological texture (antitexture) in this 
system, L w is the characteristic texture size. During the 
phase ordering process, because isolated textures expand 
[9], textures unwind by annihilating with adjacent an- 
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titextures at the boundaries between regions of positive 
and negative winding 

We begin with a heuristic argument relating L and L w . 
Consider a region of length I. If the phase angles in the 
initial condition have only short-range correlations, with 
correlation length £ , then the initial net winding over the 
length I is of order (l/^o) 1 ^ 2 - Because the total winding 
is a topological invariant, the net winding on scales much 
larger than the phase coherence length, I » L, will be 
unchanged. At later times, the length I contains of order 
l/L sections, each winding in a given sense, with of order 
L /L w windings per section. The net winding is therefore 
G/£o) 1/2 ~ {L/L w )(l/Lf/\ giving L w ~ (££ ) 1/2 . (We 
assume here that the fluctuations in the total winding per 
section are comparable with the mean winding). This 
indicates that two time-dependent lengths characterize 
the system. 

We now present an explicit calculation that verifies our 
heuristic argument and gives the growth laws for L and 
L w . It is convenient to formulate the problem in terms of 
the U(l) model for a complex scalar field <p = pexp(i6). 
We take the conventional Ginzburg-Landau free-energy 
functional 

F[<t>] = J dx [{d x <t>){d x <t>*) + (g/2)(l - W*f] . (1) 
The purely dissipative equation of motion is 

(2) 

where fi — and 2 for nonconserved and conserved dy- 
namics respectively. It is mathematically convenient to 
take the limit g — > oo, which imposes the constraint 
\4>\ = 1, corresponding to a non- linear sigma model or 
"hard-spin" description. This limit is taken by writing 
4> = cxp(i9 — (3 / g) in (2), expanding in (3/g, and retain- 
ing only terms of order unity. [ This method can be quite 
generally applied to conserved vector systems by taking 
<j> = exp(-0/g)<j>, where \4>\ = 1.] For the ID XY model, 
we obtain 

i9exp{i9) = {-d 2 x Y /2 [iO" - {6' f + 2/3] cxp(ifl) , (3) 

where dots and primes indicate derivatives with respect 
to t and x respectively. 

Consider first the nonconserved case, [i = 0. Equating 
real and imaginary parts in (3) gives 

9 = 9" , (4) 
= ■ (5) 

Thus the phase equation (4) decouples from the ampli- 
tude equation (5), and the amplitude is slaved to the 
phase. 

For the conserved case, \i — 2, the same treatment 
yields coupled equations for 9 and (3. Motivated by Eq. 
(5), we put (3 = (9') 2 /2 + 7, where we anticipate that 7 



will be negligible at late times. The resulting equations 
are 

9 = (9' f9" - 9"" - 2 7 6>" - 4 7 '6»' , (6) 
2i' = 2{9') 2 1 + 29'9"' + {9"f . (7) 

It is easy to show explicitly that these equations conserve 
the order parameter, i.e. d t J dx exp(i9) = 0. 

As a consequence of the two length scales in the prob- 
lem, the first term on the right of (6) dominates at 
late times, so that the phase equation (6) again de- 
couples from the amplitude equation (7) at late times. 
The key point is that while the typical size of 9' is 
given by 9' ~ l/L Wl the spatial variation of 9' oc- 
curs on the longer scale L (see Figure 2). Thus each 
higher derivative generates an extra factor of l/L, giv- 
ing 9" ~ l/LL w , 9"' ~ l/L 2 L w , etc. Thus on the 
right of (6) 9"" ~ l/L 3 L w is negligible compared to 
(9') 2 9" - l/LL 3 w . Now look at Eq. (7). Demanding 
that (#') 2 7 - &&" - (9") 2 - 1/L 2 L 2 W gives 7 - l/L 2 
(and on the left, 7" ~ l/L 4 is negligible). Putting this 
in (6), we find that the terms involving 7 are both of or- 
der 1/L 3 L W and therefore negligible at late times. Thus 
the first term on the right of (6) dominates at late times, 
giving the simplified dynamics 

9 = (9') 2 9" . (8) 

This equation is one of the central results of the paper, 
and represents a significant simplification of the original 
equation of motion. Although equation (8) no longer 
conserves the order parameter at all times, the omitted 
terms on the right of (6) are of relative order L 2 w jL 2 <~ 
l/L, and the conservation is asymptotically recovered at 
late times. 

The key step in deriving the growth exponents for L 
and L w is to transform from the phase variable 9 to the 
phase gradient y = 9'. Note that q(x) = y(x)/2n is 
just the local winding rate or the 'topological charge den- 
sity' at point x. Our basic assumption is that, whereas 
the order parameter representation sketched in Figure 
1 can never be made scale invariant due to the two 
different length scales, the same morphology in the y- 
representation of Figure 2 is scale invariant under a si- 
multaneous rescaling of x by L and y by 1/L^ <~ 1/L X 
(where we anticipate \ = 1/2 from our heuristic argu- 
ment) . This is an important generalization of the stan- 
dard dynamical scaling hypothesis, and is confirmed by 
exact calculation for non-conserved dynamics and by sim- 
ulation for conserved dynamics [7]. 

For compactness, we combine (4) and (8) as the sin- 
gle equation 9 = (9' 2 ) tl ^ 2 9". In terms of y, this reads 
V = [(y 2 )^ 2 !/']' 7 (y = Q')- Making the scale trans- 
formations x — > bx, t —> b z t, y — > b~ x y and demanding 
scale invariant behavior, gives our main result: 

z = 2 + nx ■ (9) 
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To determine the exponent x> an d to exemplify the 
scaling, it is convenient to introduce the squared phase 
difference correlation function 

G(r,t) = ([6(x + r,t)-9(x,t)} 2 ) 

= L 2 ^g{r/L), (10) 

where the scaling form follows from the scaling trans- 
formations above and from noting that phase differences 
scale as L 1_x . Alternatively, we could work with the 
phase gradient (or topological charge density) correlation 
function H(r, t) = (y(x + r)y(x)) = L- 2 *h{r/L). 

The angled brackets in (10) represent an average over 
an ensemble of initial conditions. A natural choice of 
initial conditions is the gaussian distribution 

P[9(x, 0)] cx exp{- 0fc(O)0-fc(O)/2<7 fe } , (11) 
k 

where #fc(0) is the Fourier amplitude of 6(x, 0). The pair 
correlation function for the order parameter at t = 
then takes the form C(r,0) = (<j)(x + r, 0)<j>*(x, 0)) = 
cxp{-([0(x + r,0) - 9(x,0)] 2 )/2} = exp{-£ fe a fc (l - 
cos/cr)} . Choosing = 2/^k 2 yields C(r, 0) = 
cxp(— r/^o), appropriate to a quench from a disordered 
phase with correlation length £o- This corresponds to 
a 'random walk' of the initial phase angles, so that 
G(r,0) = 2r/£ . 

Now consider the dynamics. Since the topological 
charge is locally conserved, the development of phase co- 
herence at scale L due to texture-antitexture annihilation 
does not affect the phase-difference correlation function 
at larger separations r ^> L, i.e. G(r,t) — > 2r/£o in this 
limit. Hence from Eq. (10), the scaling function g(x) ~ x 
for x — > oo, and since L must drop out in this limit we 
have 

X=l/2. (12) 

Putting this into (9) gives z = 2 + ^/2, and so 

L~t 2 /^\ (13) 
L w ~ L 1 ' 2 ~ iV(*+f0 . (14) 

Previous studies of phase ordering systems have usu- 
ally concentrated on the pair correlation function C (r, t) , 
and its Fourier transform, the structure factor S(k,t). 
In the present context, the phase difference correlation 
function G(r, t) is more appropriate, as it reveals the 
full scaling structure (10). The scaling properties of 
C(r, t) can, however, be inferred. The usual pair cor- 
relation function is given by C(r, t) = {exp{i[9{x + r, t) — 

9(x,t)}}) = (exp{i(ry+r 2 y' /2 + r 3 y" /6H }}, where the 

second equality follows from the Taylor series expansion 
of 9{x + r,t). In the late-time limit r —> oo, L w — > oo, 
with r/L w fixed, only the leading term in the expan- 
sion survives, because ry ~ r/L w is of order unity, while 



r 2 y' ~ r 2 / LL W is of order L w /L <C 1, and the higher 
terms are smaller still. This limit probes correlations on 
the scale L w , since L" 1 sets the scale of y = 6', so that 

C(r,t) = (exp(iry))=f(r/L w ), (15) 

where L w ~ i 1 / 4 for [i = and t 1 / 6 for [i = 2. Because 
the structure factor S(k,t) is the spatial Fourier trans- 
form of C(r,t), from (15) we see that S(k,t) = P(k,t), 
where P{y,t) is the single-point probability distribution 
for y. For = 0, the linear dynamics (4) combined 
with the gaussian initial condition (11) ensures that the 
probability distribution, and hence S(k,t) and C(r,t) 
are gaussian at all times. The exact solution of the 
model [6,7] confirms this feature, with the expected scale 
length L w ~ i 1 / 4 . For the conserved case, the con- 
servation requires that S(k,t) vanish at k = 0, imply- 
ing P(0, t) = in the scaling limit. Thus P{y, t) can- 
not be gaussian in this case, but must have a double 
peaked structure. Numerical studies [7,8] are consistent 
with the result L w ~ i 1 / 6 derived above, and indicate 
[7] that P(y,t) is approximately described by the form 
P(y, t) - L 3 w y 2 exp(-const y 2 L 2 J. 

Time-derivative correlations probe the scaling proper- 
ties of the full two-time correlations C(r,t,t'). In the 
same limit r — > oo, L w — > oo, with r/L w fixed, we 
have T(r, t) = d t d t ,\ t = v C(r,t,t') = (9 2 vcp{i[9(x + r, t) - 
9(x,t)]}) = {y' 2 y 2fJ - cxp(iry)) , where we have used equa- 
tions (4) and (8) for 9. For the non-conserved case, be- 
cause the variables are gaussian and (yy 1 ) — {(y 2 )')/2 = 
0, then T(r,t) = (y /2 ) (exp(iry)) = (y' 2 )C(r,t). For the 
conserved case, the phase variables are not gaussian so 
T{r,t) is not simply proportional to C(r, i). In both 
cases, we use y ~ L^ 1 and y' ~ (LL^,) -1 and the growth 
laws of equations (13) and (14) to determine 

T(r, t) = t- 2 ^+V^+Vf(r/L w ), (16) 

which breaks dynamical scaling because the time- 
dependent amplitude is not proportional to t~ 2 [2]. Be- 
cause the phase dynamics involves spatial gradients of y, 
the second length-scale L is introduced and dynamical 
scaling is broken. 

These results can be generalized to a broader class 
of correlated initial conditions which includes a conven- 
tional scaling solution. If we take at ~ k~ a in (11), we 
obtain G(r, 0) ~ r"^ 1 , provided 1 < a < 3. The require- 
ment due to local phase conservation that this form be 
recovered from the general scaling form (10) when r ^> L 
fixes 

X=(3-a)/2 ; z = 2 + M (3 - a)/2 , (17) 

where we have applied Eq. (9). For fi = we still 
have z = 2 with L - t 1 / 2 and L w - t( 3 -")/ 4 , but for 
(i = 2 we obtain z = 5 — a, giving L ~ t 1 /( 5_a ) and 
L w ~ £,(3-a)/2(5-a)_ For a = 1, we obtain G(r, 0) - lnr, 



3 



implying a power-law decay of C(r, 0). Simple scaling is 
recovered in the limit a — > 1 since x — > 1 implies that 
Lu, and L both grow in the same way with characteris- 
tic scale L ~ i 1 / 2 for [i — 0, and i ~ t 1 / 4 for /i = 2. 
These growth laws are just what we expect for a one- 
dimensional 0(2) system with simple scaling [2]. 

The essence of the 'energy scaling' approach that de- 
termines the growth laws for single-length scaling [2] 
can be used for an alternative derivation of the cen- 
tral result (9). In the hard-spin limit, the free-energy 
functional (1) becomes F[9] = J dx9' 2 . The en- 
ergy density is therefore e = {9' 2 ) ~ L~ 2 ~ L~ 2x . 
This gives the energy density dissipation rate as e ~ 
-LL-^ +1 \ However, e may be independently es- 
timated via e = ((5F/S9)9) ~ -(9" 2 (9' 2 y/ 2 ) ~ 
-L~ (2+ll) L- 2 ~ -L-P+x( 2 +f)]. Equating these two es- 
timates gives L ~ iV( 2 +Mx). 

To summarize, we have shown that the ordering dy- 
namics of the 0(2) model in d — 1 involves two charac- 
teristic length scales: L w ~ tV( 4 +M) which acts as the 
scaling length for the order parameter correlation func- 
tions, though T(r,t), and by implication C(r,t,t'), do 
not satisfy standard dynamical scaling, and L ~ t 2 /( 4 +^) 
which is the scaling length for correlations of the phase 
difference (or phase gradient). Working with the phase 
gradient, which is proportional to the topological charge 
density for this system, is necessary to provide a uni- 
fying framework. It is only by considering correlations 
of the phase differences that a simpler scaling descrip- 
tion emerges. It will be interesting to see to what ex- 
tent scaling violations, and/or a simplified description in 
terms of the topological charge density, occur in higher- 
dimensional texture systems [5]. 
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FIG. 1. A section of a system from a simulation us- 
ing non-conserved dynamics. Distance along the system is 
shown by the scale and the unit-magnitude order parameter 
is shown in the orthogonal plane. The windings of ±2tt (tex- 
tures/anti-textures) of scale L w and the clusters of monotonic 
winding of a larger scale L are evident. 




FIG. 2. The phase gradient y = 8' vs. distance is shown 
corresponding to the section of system shown in Figure 1. 
The characteristic length scales are indicated: L is the average 
distance between zeros of 9' , while 2n/L m is the characteristic 
magnitude of 9' . 
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